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Abstract —Network tomography means to estimate internal 
link states from end-to-end path measurements. In conventional 
network tomography, to make packets transmissively penetrate a 
network, a cooperation between transmitter and receiver nodes is 
required, which are located at different places in the network. In 
this paper, we propose a reflective network tomography , which can 
totally avoid such a cooperation, since a single transceiver node 
transmits packets and receives them after traversing back from 
the network. Furthermore, we are interested in identification of 
a limited number of bottleneck links, so we naturally introduce 
compressed sensing technique into it. Allowing two kinds of 
paths such as (fully) loopy path and folded path, we propose 
a computationally-efficient algorithm for constructing reflective 
paths for a given network. In the performance evaluation by 
computer simulation, we confirm the effectiveness of the proposed 
reflective network tomography scheme. 

I. Introduction 

Tomography refers to the cross-sectional imaging of an 
object from either transmission or reflection data collected 
by illuminating the object from many different directions m. 
When the object is an information network, it is called network 
tomography Gl, which has been used to encompass a class 
of approaches to infer the internal link states from end-to- 
end path measurements a. The end-to-end path behaviors 
have been transmissively measured via a cooperation between 
transmitter and receiver nodes, which are located at different 
places in a network. However if it is possible to eliminate such 
a cooperation, network tomography would become a more 
powerful method with special properties ( implementability , 
adaptability and asynchronism) for measuring and analyzing 
network specific characteristics. 

In this paper, according to the types of end-to-end path 
measurements acquisition, we first classify network tomog¬ 
raphy into transmissive and reflective network tomography, 
and after discussing their characteristics, we propose a new 
reflective network tomography scheme. Here, in the reflective 
network tomography scheme, we focus only on identification 
of a limited number of links with large delays in a network, 
where such links are referred to as bottleneck links. In this 
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Fig. 1: Transmissive end-to-end path measurement. 


scheme, a node acts as both a transmitter and a receiver, i.e., 
as a transceiver: it transmits multiple packets over a network 
along pre-determined different paths and receives the packets 
after they traverse back from the network. On the other hand, 
network tomography is formulated as an undetermined linear 
inverse problem and it cannot be always solved. However, 
the assumption in the bottleneck link identification makes it 
possible to use compressed sensing technique. To propose the 
new reflective network tomography scheme, we tackle two 
problems: how to formulate the tomography scheme and how 
to determine going around paths from/to a transceiver node. 

Note that, although end-to-end path measurements can 
be conducted either actively or passively, reflective network 
tomography scheme is only based on active measurements. 
Thus, we particularly consider active tomographic scheme in 
this paper. 


II. Network Tomography 
A. Transmissive Network Tomography 

In this subsection we define transmissive network tomog¬ 
raphy via some examples 0-||9l which are characterized by 
transmissive end-to-end path measurements. Fig. \T\ shows an 
example of a transmissive end-to-end path measurement 0. 

In a network with a defined boundary, it is assumed that 
access is available to nodes at the boundary, but not to 
any in the interior. In order to get transmissive end-to-end 


path measurements, some boundary nodes are selected as 
transmitter and receiver nodes. For example, in 0, two 
nodes are respectively assigned as a transmitter and a receiver, 
whereas in 0, there are many transmitter and receiver nodes. 
The transmitter nodes send probe packets to all (or a subset 
of) the receiver nodes to measure packet attributes on the 
paths between them. Accordingly, each probe packet trans- 
missively penetrates the network along a measurement path , 
and brings a transmissive end-to-end path measurement. In 
0 , a transmissive tomographic methodology based on unicast 
communication is proposed. In |7) and 0, on the other hand, 
a single-source multicast transmission by a single or multiple 
transmitter nodes is applied to networks with tree and general 
topologies, respectively. From such transmissive end-to-end 
path measurements between transmitter and receiver nodes, the 
internal network states such as link-level network parameters 
can be estimated. For example, in 0, link delay variance is 
estimated from transmissive end-to-end path measurements in 
a multicast setting. 

B. Reflective Network Tomography 

Unlike transmissive network tomography, reflective network 
tomography eliminates the need for special-purpose coopera¬ 
tion from receiver nodes. Namely, an end-to-end path measure¬ 
ment is calculated from records on only one node. A boundary 
node is selected as a transceiver node, and it injects probe 
packets into the network. Each probe packet goes and back to 
the transceiver node along a different measurement path, and 
brings reflective end-to-end path measurements. For example, 
in da, a reflective network tomography scheme based on 
round trip time (RTT) measurement only along a folded path 
(see IIV-BI for its definition) is proposed to estimate the delay 
variance for a link of interest. Thus, in contrast to transmissive 
network tomography, reflective network tomography is defined 
by reflective end-to-end path measurements. 

III. Properties of Reflective Network 
Tomography 

A. Implementability 

The methods described in the above transmissive network 
tomography all require a coordination between transmitter 
and receiver nodes. However, the following problems have 
not been discussed deeply: how to access all the transmitter 
and receiver nodes and how to establish the coordination 
between them, in order to implement the network tomography, 
i.e., designate the measurement paths, transmit active probe 
packets and collect the end-to-end path measurements. In a 
network, these would occupy some part of the time/frequency 
resource and consume some energy. Without any solution 
strategy, these problems would limit the scope of the paths 
over which the measurements can be made. Thus most of 
them would not be widely applicable because of the lack of an 
available widespread infrastructure for transmissive end-to-end 
path measurements. 

On the other hand, the reflective network tomography 
scheme does not require special cooperation from the other 


interior and boundary nodes, because the reflective end-to-end 
path measurements are calculated only by a single transceiver 
node. We just use the transceiver node to implement the re¬ 
flective network tomography, so we can say that the reflective 
network tomography can be carried out more easily. 

B. Adaptability 

Most of the existing transmissive network tomography 
schemes are based on non-adaptive measurements in them¬ 
selves. Namely, the measurement paths are often fixed in 
advance and do not depend on the previously acquired mea¬ 
surements. The reason is that it is difficult to feed back the 
prior end-to-end path measurements from receiver nodes to 
transmitter nodes every probing. 

In the reflective network tomography scheme, on the other 
hand, since the probe packets return to the transceiver node, 
measurement paths can be adaptively selected depending on 
the previously gathered information. So it can give us the 
advantage of sequential measuring schemes that adapt to 
network states using information gathered throughout a mea¬ 
surement period. Furthermore, many current methodologies 
usually assume that network states are stationary throughout 
the tomography period. Even when this assumption is not 
satisfied, however, reflective network tomography scheme may 
be workable thanks to its adaptability. 

C. Asynchronism 

When focusing on transmissive delay tomography which 
is transmissive network tomography for link delays, end- 
to-end path measurements are usually calculated from the 
transmission time and reception time reported by the trans¬ 
mitter and receiver nodes, respectively. Therefore, it requires 
clock synchronization between them. However, the clock 
synchronization is sometimes hard to achieve or not guaran¬ 
teed, especially in wireless networks such as wireless sensor 
networks, in which electronic components of nodes are too 
untrustable to meet the requirement of clock synchronization 
in terms of accuracy and complexity (TO, G2. So, although 
delay tomography scheme workable in clock-asynchronous 
networks is preferable, to the best of the authors’ knowledge, 
the transmissive synchronization-free network tomography has 
been studied only in (H. 

On the other hand, reflective network tomography scheme 
does not require any clock synchronization for any other 
nodes in a network. The time delay for a packet traveling 
through a measurement path can be estimated by checking the 
transmission time and reception time on a transceiver node’s 
clock. Therefore, the reflective network tomography scheme 
is potentially available in clock-asynchronous networks. 

IV. Proposed Reflective Network Tomography 
Scheme 

A. Compressed Sensing 

Compressed sensing is an effective theory in signal/image 
processing for reconstructing a finite-dimensional sparse vec¬ 
tor based on its linear measurements of dimension smaller than 




the size of the unknown sparse vector lfl4l . |[T5ll . Recently, 
compressed sensing has been also used for network tomogra¬ 
phy 0|, 0, lff6l . In this subsection, as the preliminary for 
compressed sensing, we give several definitions. 

First, we define the i p norm (p > 1) of a vector x = 
[x\ X 2 • • • xj] T G 7 Z J as 

iwip = (S w j P > (1) 

i=l 

where T denotes the transpose operator. 

Next, we assume that, through a matrix A G 7 Z lxj (I < J), 
we obtain a linear measurement vector y = [yi 1/2 ■ • ■ yi] T G 
VJ for a vector x = [x\ X 2 • • • xj] T G 1Z J as y = Ax. 
Whether or not one can recover a sparse vector x from y by 
means of compressed sensing can be evaluated by the mutual 
coherence /i(A) lfl5l . To calculate the mutual coherence of 
A, by picking up the j'-th and j’-th column vectors from A 
we construct the partial matrix as 

A jj' = [c j C y], (2) 

where c j and Cj> are the j-th and j'-th column vectors of 
A, respectively. The mutual coherence /i(A) is defined as the 
maximum value of v(Ajj') (1 <j,f < jj 7 
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then there exists at most one vector x with at most k nonzero 
components that y = Ax. 

B. System Model 

We consider a delay tomographic scheme which identifies 
a few bottleneck links in an asynchronous network from re¬ 
flective end-to-end path measurements. Our approach employs 
unicast communication. Let Q = (V,£) denote an undirected 
network, where V is the node set, and £ C V x V is the link 
set. Note that (i, j) G £ implies that (j, i) G £ since the graph 
is undirected. We assume that the topology is fixed throughout 
the measurement period and there is only one transceiver node 
5. 

Due to the fact that the overall delay of a path is the sum of 
the delays of all links belonging to the path, delay tomography 
problem can be formulated as an inverse problem to recover 
link delays based on linear measurements. Here, we measure 
packet traveling times (PTTs) along two kinds of paths by 
injecting probe packets into the network. One is a (fully) loopy 
path (LP) defined as the one where any nodes do not appear 
more than once except for a transceiver node, and the other is 
a folded path (FP) defined as the one where any nodes appear 

Actually our proposed scheme can also be extended to directed graph 
models. 



twice except for a destination node (in other words, an FP is 
composed of a path from a transceiver node to a destination 
node and a path from the destination node to the transceiver 
node along the same undirected curve between them). Fig. [2] 
shows an example for LP and FP, and we do not consider 
any path containing partial loops ( routing constraint). Now, 
we define W = {path^\ l = 1 , 2 ,..., |W|} as a subset 
of all paths to measure PTTs based on s , where path ® = 
{(s, (v( l,1 ),v( 1,2 )),..., (t/ , ’IP a ^ i) l- 1 ), s)} C £ repre¬ 
sents the Z-th path in W and G V \ {s} (m = 

1,..., | path^p | — 1) are intermediate nodes in the path. 

We reformulate W and £ as W = {wi,W 2 ,... ,wi} and 
£ = {ei, e 2 ,..., ej}, respectively, where / = |W| and 
J — \£\ denote the numbers of paths and links, respectively. 
We assume that link delays d ej arise independently on each 
link ej (j = 1, 2,..., J), which does not depend on the 
direction. Thus, a probe packet transmitted on a path Wi 
(i = 1, 2,...,/) is successfully returned to s with total delay 
D w . = XL .a,,, d e <- We define measurement vector y = 
[yi y 2 • • • yi] T and link delay vector x = [x\ X 2 • • • xj] T as 


yi — D Wi — ^ ^ d e ., 

ejEwi ( 6 ) 

Xj = d e .. 

Then, we obtain 

y = Ax, (7) 


where A G {0,l,2} /xJ represents the (reflective ) routing 
matrix of W, i.e., (i, j)-th component (i = 1, 2 ,..., /, j = 
1, 2,..., J) in A is set to = 1 or — 2 if ej G Wi , and 

aij = 0 otherwise. The row size I is related to the interval 
devoted to the tomography scheme, and the entrywise matrix 
norm of A 

! a II = (EEM (8) 

i= 1 3 = 1 

is related to the traffic load of probe packets. The I and || A|| 
determine the energy required for accomplishing a tomography 
scheme, so the former is referred to as the interval factor , 
whereas the latter the traffic factor. For a given detectability 
of bottleneck links, the two factors of a better routing matrix 
should be smaller. 

Note that link states are assumed to be stationary, i.e., link 
delays do not change while the proposed scheme is applied, 




TABLE I: Symbols in Algorithm [T] 


P disjoint 

Set of node-disjoint paths. 

P disjoint 

Set of node-disjoint reverse paths, which are constructed by reversing directions of paths in Pdisioint 

w 

Set of definitive measurement paths. 

path ^ + path ^ 

Path connecting path ^ and path^ b h 

Pall 

Set of all candidates for measurement paths. 

A) 

Function which returns the mutual coherence of A if no column vector equals 0, and number greater than 1 otherwise. 

getCostMin( / P) 

Function which returns a path whose cost is the minimum in a path set V. 


and a few bottleneck links exist in the network. Next, because 
it is possible to approximate the elements of x corresponding 
to small link delays to be zero by attributing the delays only to 
the few bottleneck links, the idea of compressed sensing can 
be naturally introduced to network tomography. So we utilize 
compressed sensing based on optimization [17], [18] in 
order to reduce traffic load of probe packets. Finally, when 
using the PTTs, the assumption of the undirected graph may 
lead to inaccurate estimates given the asymmetric communi¬ 
cation HE). However, we are interested in identification of a 
limited number of bottleneck links, thus, the assumption can be 
considered to be valid, since it does not require measurements 
with accuracy. 

C. Routing Matrix Construction 

Now, we propose a simple algorithm composed of two 
steps for constructing a routing matrix A. This algorithm is 
for a reflective routing matrix which can quickly identify a 
bottleneck, assuming that a bottleneck link rarely arise in 
the network. Algorithm [T] shows the algorithm, and Table U 
describes symbols used in Algorithm [T] 

First, in STEP 1 the algorithm constructs a set of 
paths as measurement path candidates based on node- 
disjoint paths algorithm described in [20]. The function 
NodeDisjomtAlgorithm(,s,?;) in Algorithm [I] returns the maxi¬ 
mum set of node-disjoint paths from s to r. The set of node- 
disjoint paths implies the shortest combination of paths where 
no nodes are shared among the paths. By connecting every 
node-disjoint path from 5 to v (for all v G V\s), this algorithm 
lists up the candidates for measurement paths, which satisfy 
the routing constraint that any path is an LP or an FP. 

Then, out of the path candidates constructed by STEP 
1, STEP 2 selects paths as measurement paths one-by-one 
according to the cost of candidates until the mutual coherence 
of the constructed routing matrix becomes less than 1.0. If 
several paths have the same minimum cost, the shortest path 
is selected out of them. Here, we define the cost function for 
a measurement path (path G Pail) as 

Cost {path) 

(Number of unused links in W out of path ) -1 

(if -f)i(A) > 1.0) 

Number of = 1 (1 < j,j' < J:j ^ /) 

(otherwise), 

where A' is constructed from a set W + {path}, and this cost 


function is used in getCostMin(P). Once the mutual coherence 
of the constructed matrix becomes less than 1.0, this algorithm 
terminates. STEP 2 cannot directly select a path depending on 
the number of nodes over the path. Therefore, the proposed 
routing matrix construction algorithm pays attention to the 
interval factor rather than the traffic factor. 

Algorithm 1 Proposed Routing Matrix Construction Algo¬ 
rithm_ 

Require: Network Topology and s. 

Ensure: Routing Matrix A. 

STEP 1 : Search for path candidates 
for all v G V \ s do 

Pdisjoint •= NodeDisjointAlgorithm^, v ). 
for all path ^ ^ G Pdisjoint(p — 1, 2 • • • |Pdisjoint|) do 
for all path m e Pdisjoint(& “ 1 ,2--- |Pdisjoint|) do 
Vail '■= Van U {path^ + path (b) }. 

end for 
end for 
end for 

STEP 2 : Selection of measurement paths 
while T)x(A) >1.0 do 

path^ min ^ := getCostMin (V a \\ \ W). 

W:=>VU {path (min) }. 

Construct A from set W. 
end while 
return A. 


V. Performance Evaluation 

In this section, we discuss the following items: 

• Can the proposed algorithm (Algorithm |T]) construct a 
fully adequate routing matrix? 

• Can a routing matrix constructed by the proposed al¬ 
gorithm actually identify a bottleneck link in a network 
where only one bottleneck link exists? 

• How does a routing matrix with smaller interval and traf¬ 
fic factors behave in a network where several bottleneck 
links exist? 

Fig. Eta) shows the network topology with 8 nodes and 11 
links used for performance evaluation by computer simulation, 
where there is only one transceiver node s. We assume that the 
delay of a bottleneck link is constant with x^ whereas that of 
a normal link denoted by x( N ) is independent and identically 
distributed (i.i.d.) with average a x m and standard deviation 
cr x (N). In this paper, we also assume that all the nodes are 
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Fig. 3: Network topology (8 nodes and 16 links) and the measurement paths constructed by the proposed algorithm. 


TABLE II: Routing Matrices 


Matrix 

Size 

Entry wise 
Matrix Norm 

LP:FP 

Numbers 

Number of 
Candidates 

Pi 

5 x 11 

30 

3 : 2 

18paths 

P 2 

5 x 11 

28 

4 : 1 

18paths 

Pa 

4 x 11 

28 

3 : 1 

78paths 


wirelessly connected so is Gaussian-distributed [21] with 
a x (N) = 15 msec and a x ( N ) = 3 msec El, El. 

First, Table [III shows the constructed three routing matrices 
whose mutual coherences are less than 1. In Table HU Pi is 
constructed by the proposed algorithm composed of STEP 1 
and STEP 2 (the measurement paths are shown in Figs. Hb)- 
(f)), P2 is constructed by a greedy search from the path 
candidates listed by STEP 1 (instead of STEP 2, the paths are 
selected from all combinations of the path candidates by STEP 
1 , which minimizes the interval and traffic factors), and P3 is 
also constructed by a greedy search from path candidates listed 
by STEP 1 and additional FP candidates (all FPs are added to 
the path candidates by STEP 1 and then the paths are selected 
from all combinations of the increased path candidates, which 
minimize the interval and traffic factors). It is impossible for 
the proposed algorithm to always select the paths which really 
minimize the interval and traffic factors due to its one-by-one 
policy (in STEP 2), on the other hand, the greedy search- 
based algorithms can always select the optimum set of paths 
from all combinations of paths. Comparing Pi and P2, the 
proposed algorithm composed STEP 1 and STEP 2 constructs 
the routing matrix whose traffic factor is a little larger, and 
comparing P 2 and P 3 , the number of path candidates by STEP 
1 seems insufficient. However, when the size of network is 
large, for the case where I measurement paths are selected 
from N candidates, the proposed algorithm calculates the 



Fig. 4: Bottleneck link detection ratio vs. bottleneck link delay 
X ( B ) (Number of bottleneck links k = 1). 


cost function (IN + 1(1 — l)/ 2 ) times, whereas the greedy 
search-based algorithms lead to combinatorial explosion. So, 
taking into consideration that computational complexity of the 
proposed algorithm is much lower than that of the greedy 
search-based algorithm, it can be concluded that the proposed 
algorithm can efficiently construct a fully adequate routing 
matrix. 

The termination of the proposed algorithm is guaranteed 
since as the number of measurement paths increases, the mu¬ 
tual coherence of the constructed routing matrix monotonously 
decreases. While mutual coherence can provide a guarantee 
of the recovery of exactly sparse vectors, the link delay 
vector x is approximately sparse in the model for performance 
evaluation. Therefore, to confirm whether or not the bottleneck 
link detectability of the reflective network tomography scheme 
is consistent with the meaning of the mutual coherence. So, we 
assumed that there is a bottleneck link in the network, that is, 
we set the number of bottleneck links k to 1 in the computer 






































Fig. 5: Bottleneck link detection ratio vs. number of bottleneck 
links k (Bottleneck link delay X (B) = 1.0). 


Some technical issues remain in the proposed scheme. 
First, we have to propose a better routing matrix construction 
algorithm, and evaluate reflective network tomography on 
larger networks. We also have to propose an adaptive network 
tomography to take advantage of reflective characteristic. 
Since these issues are beyond the scope of this paper, we leave 
them as future works. 
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simulation. Here, we also define bottleneck link detection ratio 
which is defined as the number of correctly detected bottleneck 
links divided by the total number of given bottleneck links. 
Fig. [4] shows the bottleneck link detection ratio versus the 
bottleneck link delay for k = 1. Although the link delay 
vector x is not exactly sparse, as the bottleneck link delay 
x becomes larger, the bottleneck link detection ratios of 
the three routing matrices approaches 1.0. This means that, if 
the bottleneck link delay x^ is fully larger, the link delay 
vector x can be regarded approximately as a sparse vector, 
and the mutual coherence can also guarantee the recovery of 
approximately 1-sparse vectors. Thus, it can be concluded that 
the proposed scheme can effectively detect a bottleneck link. 

Finally, we evaluated the matrices Pi,P 2 ,P 3 in the net¬ 
work with the number of the bottleneck links k = 1,2,3. 
Fig. 0 shows the bottleneck link detection ratio versus the 
number of bottleneck links, where we set x^ to 1 sec, 
which corresponds to about 66.6 times as large as ol x {n). 
For k > 2, all the bottleneck link detection ratios fall down 
sharply. This is because that the algorithms introduced here 
all try to construct routing matrices with smaller interval and 
traffic factors, which have worse impact on the bottleneck link 
detectability of k > 2 (the proposed algorithm composed 
of STEP 1 and STEP 2 pays attention only to reducing 
the interval factor, but it also results in reduction of the 
traffic factor). Therefore, for a network with the possibility 
that several bottleneck links arise simultaneously, we need to 
redesign the termination condition and cost function. 

VI. Conclusion 

In this paper, according to the types of end-to-end path 
measurements acquisition, we classified network tomography 
into transmissive and reflective schemes and proposed a new 
reflective network tomography with their advantageous charac¬ 
teristics over conventional transmissive network tomography. 
We proposed a simple reflective routing matrix construction 
algorithms composed of two steps, and by computer simulation 
we showed that it can effectively construct an adequate routing 
matrix guaranteeing a designed bottleneck link detectability of 
km 1 . 
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